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Abstract 


Temporal and spatial variation of fine-structure constant a = e 1 2 /he in cosmology has 
been reported in analysis of combination Keck and VLT data. This paper studies 
this variation based on consideration of basic spacetime symmetry in physics. Both 
laboratory ao and distant a z are deduced from relativistic spectrum equations of atoms 
(e.g.,hydrogen atom) defined in inertial reference system. When Einstein’s A / 0, the 
metric of local inertial reference systems in SM of cosmology is Beltrami metric instead 
of Minkowski, and the basic spacetime symmetry has to be de Sitter (dS) group. The 
corresponding special relativity (SR) is dS-SR. A model based on dS-SR is suggested. 
Comparing the predictions on a-varying with the data, the parameters are determined. 
The best-fit dipole mode in a’s spatial varying is reproduced by this dS-SR model, a- 
varyings in whole sky is also studied. The results are generally in agreement with the 
estimations of observations. The main conclusion is that the phenomenon of a-varying 
cosmologically with dipole mode dominating is due to the de Sitter (or anti de Sitter) 
spacetime symmetry with a Minkowski point in an extended special relativity called de 
Sitter invariant special relativity (dS-SR) developed by Dirac-Inonii-Wigner-Giirsey- 
Lee-Lu-Zou-Guo. 
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1 Introduction 


Temporal and spatial variation of fundamental constants is a possibility, or even a neces¬ 
sity, in an expanding Universe (see, [1], and a review of 0). A change in the fine structure 
constant a = e 2 /(hc) could be detected via shifts in the frequencies of atomic transitions in 
quasar absorption systems. Recent analysis of a combined sample of quasar absorption line 
spectra obtained using UVES (the Ultraviolet and Visual Echelle Spectrograph) on VLT (the 
Very Large Telescope ) and HIRES (the High Resolution Echelle Spectrometer) on the Keck 
Telescope have provided hints of a spatial variation of the fine structure constant (a), which 
is well represented by an angular dipole model 0 ffl. That is, a could be smaller in one 
direction in the sky yet larger in the opposite direction at the time of absorption. Prior to 
that, measurements of possible time variations of the fine-structure constant was achieved by 
the same method in the Keck telescope. It has been shown that the time variation of a 
exists. Thus, a is a “constant” varying with both red-shift z (or cosmologic time) and direc¬ 
tion in the sky equatorial coordinates. Namely, a = a z (fl) where fi indicates the direction 
in the sky. Those direct measurements of possible space-time variations of the fine-structure 
constant are of utmost importance for a complete understanding of fundamental physics. 

A straightforward conjecture for this phenomenon is that the space-time function of a(x) 
may be thought as a scalar field (p{x ) or a function of (p{x ) in the spacetime with some 
suitable dynamics (see, e.g, m mi na ). The ip(x) is a matter field and fills the Universe 
everywhere. Sometime one could call it dilaton-like scalar field. Along this way of thinking, 
authors of reference [13] argued that the spatial variation of the fine structure constant a 
may be attributable to the domain wall of (p(x) in the Universe. 

In this present paper, we would like to present a matter-field-free scheme to answer the 
challenging questions such as why the fine-structure “constant” a varies over space-time, 
and why spatial variation of a is well represented by an angular dipole mode. The scheme 
is still in the framework of standard cosmology and of the Special Relativity (SR) theory 
except that SR’s spacetime symmetry will be extended. Concretely, we shall apply the de 
Sitter invariant Dirac equation to the distant hydrogen atom to explain such variations of 
a in cosmology. The calculations are based on the theory of the de Sitter invariant special 
relativity (dS-SR) developed by Dirac-Inonu-Wigner-Gursey-Lee-Lu-Zou-Guo [141122] . To 
study atom physics in dS-SR were firstly called for by P.A.M. Dirac in 1935 [14] . 

To show the fine-structure constant a is unvarying over space-time in the Standard Model 
(SM) of physics including cosmology, we examine the relativistic wave equation of an electron 
in hydrogen in SM of physics. First, we consider the laboratory atom. From the viewpoint 
of cosmology, the energy level E of a free hydrogen atom in laboratory is determined by the 
Dirac equation in a local inertial coordinates system located at the Earth in the Universe 
described by Friedmann-Robertson-Walker (FRW) metric. The spacetime metric of the local 
inertial system is Minkowski metric: 
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which is spacetime independent. E satisfies Dirac spectrum equation: 


Eif> = \^—ihcot • V —-—b m e c 2 /3J ip, (1.2) 

where {cc = a 1 i + a 2 j + a 3 k, /3} are Dirac matrices, V = (d/dx 1 )i+(d/dx 2 )j + (d/dx 3 )k and 
r = yj (x 1 ) 2 + (x 2 ) 2 + (a; 3 ) 2 . This matrix-differential equation is integrable and the solution 
of the eigenvalue E is (see, e.g., |23] ) 

- 1/2 

(1.3) 


E = W n « = m e c 1 + 
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hr 
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k | + s) 2 / 

«| = (j + 1/2) = 1, 2, 3 
n — 1, 2, 3 • • • . 


We keep in mind that the coefficients of operators — icn ■ V and — 1/r in Eq. (ll.2p are he and 
e 2 respectively, and their ratio is the definition of a (see Eq (jl .3(1 1. 

Next, we consider a distant atom of hydrogen located on the light-cone of FRW-Universe 
(see FigH]), i.e., the nucleus coordinate is Q M (z) = {Q°, Q}, and electron’s is L^(z) = 
{L°, L}. Noting that the metric of the local inertial coordinate system at Q M in FRW- 
Universe is still r] fll , (II.ip because of the spacetime-independency of r] fll , and denoting L^(z) — 
Q^(z) = x ,fl , the electron wave equation in the distant atom reads 


E'ip = ( —ihcoL ■ V 


' - e ~ + m e c 2 j^j ip, 


(1.4) 


where V' = (d/dx n )i + ( d/dx' 2 )} + ( d/dx ,3 )k and r' = yf (x' 1 ) 2 + (a;' 2 ) 2 + (x' 3 ) 2 . Then we 
find out that 

„2 


a = a 7 = 


he 


Comparing (11.51) with (II,3p . we conclude that 


a, = ct, 


(1.5) 


( 1 . 6 ) 


which indicates that the fine-structure constant is unvarying indeed in SM of physics. 

The argument above on a-unvarying in SM made by comparing Dirac equation for lab¬ 
oratory hydrogen atom with that for a distant hydrogen atom is deeply related to aspects 
of Special Relativity (SR), General Relativity (GR) and Cosmology. In other words, the a- 
varying phenomena reported in [3H8] implies some new physics beyond SM. Further remarks 
on this issue are follows: 

1. The spectrum equations (1 1.2 p (II.4p come from the following inhomogeneous-Lorentz 
(or Poincare) invariant (i.e., ISO( 3,1)) Dirac equation and Maxwell equation in local 
inertial systems of FRW Universe: 

,L m eC, 


iyfD 1 . 


1 )ip = 0, 


(1.7) 

( 1 . 8 ) 


" h 

F»\ v = f = -(P 0 47recf (3) (x), 

where D^ —ie / (ch)rj flu A u , and the electromagnetic potential A u = {0 = e/r, A}. 

So, the operator structure of (1 1.2 j) . (11.41) and a dimensionless combination of universal 
constants a = e 2 /{he) are rooted in the symmetry assumption of the theory. 
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2. The point that a is unvarying is deduced from the constancy of the adopted metric 
of local inertial coordinate system in the FRW Universe (i,e., {rj^} =const.). So, the 
fact of the a-varying in real world reported in [3HE] indicates that the metric of local 
inertial coordinate system in the real Universe may be spacetime-dependent. 

3. Minkowski metric rj lw = diag{+, —, —, —} is the basic spacetime metric of Einstein’s 
Special Relativity (E-SR) in SM. The most general transformation to preserve metric 
rj pu is Poincare group (or inhomogeneous Lorentz group ISO(l,3)). It is well known 
that the Poincare group is the limit of the de Sitter group with pseudo-sphere radius 
\R\ — > oo. Therefore, E-SR may possibly be extended to a SR theory with de Sitter 
space-time symmetry. Since P.A.M. Dirac’s work in 1935 031 many discussions (e.g., 
E. Inonii and E. P. Wigner in 1968 |T5]; F. Giirsey and T.D. Lee in 1963 |T6], etc ) 
pointed to such a possible extension of E-SR. In 1970’s, K.H.Look (Qi-Keng Lu) and 
his collaborators Z.L.Zou, H.Y.Guo suggested the de Sitter Invariant Special Relativity 
(dS-SR) [lTj jUj (see Appendix A, and also [T9I120] and Appendix in [M] for the English 
version). It has been proved that Lu-Zou-Guo’s dS-SR is a satisfying and self-consistent 
special relativity theory. In 2005, one of us (MLY) and Xiao, Huang, Li suggested dS- 
SR Quantum Mechanics (QM) [20 ]. 

4. Beltrami metric (see Appendix A) 

Bn V (x) = rjfjv/a(x) + 'n^x x r] up x p /(R 2 a(x) 2 ), with a(x) = 1 - r) p „x p x u /R 2 > 0 (1.9) 


is the basic metric of dS-SR with Minkowski point coordinates M p = 0 (i.e., B fiu (x )\ X=M=0 
= Vnv) EHJ. Both and B pu lead to the inertial motion law for free particles x = 0, 
which is the precondition to define inertial reference systems required by special rel¬ 
ativity theories. However, the LQ,,-preserving coordinate transformation group is de 
Sitter group SO (A, 1) (or £0(3, 2)) rather than E-SR’s inhomogeneous Lorentz group 
ISO( 1,3) jT7H2U] . different from the case of rj lw . In addition, r] lw does not satisfy the 
Einstein equation with A (Einstein cosmology constant) in vacuum, but B pv (x) does 
satisfy it (see below). Generally, when ^ 0, the basic metric of dS-SR is modified 
to be 


BW(x) = B^(x - M) 


rip, %x(x^ ~ - M") 

cr( M \x) R 2 <j( M ')(x) 2 


where 


a^ M \x) = a(x — M) 


Vuvix 11 — M p )(x v — M u ) 
R 2 


( 1 . 10 ) 

( 1 . 11 ) 


which will be called Modified Beltrami metric, or M-Beltrami metric. Based on 
B^\x), the dS-invariant special relativity with M p ^ 0 can be built. The procedures 
and formulation are similar to ordinary dS-SR in [171 - 120] . which is actually a slight 
extension of usual dS-SR (see Appendix B). It is essential, however, that B$\x) is 
spacetime dependent and has more parameters { R , M p }, which may provide a possible 
clue to solve the puzzle of a-varying. 


5. Our strategy in this present paper for solving this puzzle is to pursue the following 
dS-SR Dirac equation for both electron in laboratory hydrogen and electron in distant 
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hydrogen in the FRW Universe (see Eq.(25) in [[21]): 


Til C 

- -tw = ( L12 ) 

where — ^uj a /a a b — ie/(ch)B^A v , ej 1 is the tetrad, u ab ^ is spin-connection, 

and the electromagnetic potential A v = {0 B , A}. Unlike E-SR Dirac equation (II.7j) , 
the spacetime symmetry of (II . 12(1 is de Sitter invariant group 50(4,1) (or 50(3,2)) 
instead of former 750(3,1). Specifically, the dS-SR Dirac spectrum equation can 
be deduced from (11.12j) . It is essential that the result will be different form E-SR 
equation (11.2)1 . Following the method used in (11.21) (11.31) . the coefficients of resulting 
(—icx ■ V)-type and (—l/r)-type operator terms in the dS-SR Dirac spectrum equation 
are of h z (fl)c and e z (f2) 2 respectively. Then their ratio yields prediction of a z (ft) = 
e z (f i) 2 /[h z {yi)c). The adjustable parameters in this model are R and the position of 
Minkowski point MU For simplicity, we take = {M°, M 1 , 0, 0}. It turns out to 
be a good a choice for solution to the puzzle of a-varying. 

6. Different from Quantum Mechanics (QM) wave equation (11.71) deduced from r]^, the 

equation (jl. 12j) is actually a time-dependent Hamiltonian problems in QM. This is 
because is time-dependent.Therefore the corresponding Lagrangian L d $ (see 

Eq. (lA.150|) ) and hence Hamiltonian is time-dependent [20]. In this paper, the adia¬ 
batic approach [26] [27] [28] will be used to deal with the time-dependent Hamilto¬ 
nian problems in dS-SR QM. Generally, to a H(x,t), we may express it as H(x,t) = 
H 0 ( x ) + H'(x,t). Suppose two eigenstates |s) and | m) of H 0 (x ) are not degenerate, 
i.e., A E = h(u) m — uj s ) = hoj ms ^ 0. The validity of for adiabatic approximation relies 
on the fact that the variation of the potential H'(x,t) in the the Bohr time-period 
(ATtns 0h,) )H' ms = (2 tt/ u ms )H' ms is much less than hu ms , where H' ms = (m\H'(x,t)\s). 
That makes the quantum transition from state |s) to state | m) almost impossible. 
Thus, the non-adiabatic effect corrections are small enough (or tiny) , and the adia¬ 
batic approximations are legitimate . For the wave equation of dS-SR QM of atoms 
discussed in this paper, we show that the perturbation Hamiltonian describes the time 
evolutions of the system H'(x,t) oc (c 2 t 2 /7? 2 ) (where t is the cosmic time). Since R is 
cosmologically large and R » ct , the factor ( c 2 t 2 /R 2 ) will make the time-evolution 
of the system so slow that the adiabatic approximation works. We shall provide a 
calculations to confirm this point in the paper. By this approach, we solve the station¬ 
ary dS-SR Dirac equation for one electron atom, and the spectra of the corresponding 
Hamiltonian with time-parameter are obtained. Consequently, we find out that the 
electron mass m e , the electric charge e, the Planck constant h and the foie structure 
constant a = e 1 /{he ) vary as cosmic time goes by. These are interesting consequences 
since they indicate that the time-variations of fundamental physics constants are due 
to solid known quantum evolutions of time-dependent quantum mechanics that has 
been widely discussed for a long history (e.g., see [28] and the references within). 

7. Finally, we argue that it is reasonable to assume that the Beltrami metric is the 
appropriate metric for the spacetime of the local inertial system in real world. If we 
express the total energy momentum tensor T^ u as the sum of a possible vacuum term 
—P(v)9iiv and a term T// arising from matter (including radiation), then the complete 
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(1.13) 


Einstein equation is [29] [3D] [51] : 

TZ,w - ^g^TZ + A g^u = -8nGT™ - 8vr Gp^g^, 

where p(„) is the dark energy density, so Ad ar k energy = 87 t Gp( v ), and A is originally 
introduced by Einstein in 1917, and serves as a universal constant in physics. We call 
it the Einstein (or geometry) cosmological constant. The effective cosmologic constant 
A e ff = A + Adark energy — 1.26 x 10 56 cm 2 is the observed value determined via effects 
of accelerated expansion of the universe [32] and the recent WMAP data [33]. We have 
no any a priori reason to assume the geometry cosmologic constant to be zero, so the 
vacuum Einstein equation is: 

- ^g^TZ + A g^ = 0, (1.14) 

instead of G^ = 0, and hence the vacuum solution to fll.ldp is g /t „ = with 

\R\ = \J 3/A instead of g^ = r] )W . Therefore, we conclude that the metric of the 
local inertial coordinate system in real world should be Beltrami metric rather than 
Minkowski metric. Thus, the dS-SR Dirac equation (1 1.12 p (instead of E-SR Dirac 
equation (11.71) 1 is legitimate to characterize the spectra in the real world, and then the 
a-varying over the real world space-time would occur naturally. 

This paper provides an understanding of the a-varying in cosmology reported in 0 H by 
means of extending the basic spacetime symmetry in the local inertial coordinate systems 
in the standard cosmologic model. The contents of the paper are organized as follows. In 
section 2, light-cone of Friedmann-Robertson-Walker universe is described, and relation of 
cosmological time to redshift z is shown. The relation of t — z is based on ACDM model and 
the cosmology parameters (H 0 , D m0 , Da — 1 — D m0 ) hi real world; In section 3, we describe 
the local inertial coordinate system in light-cone of FRW Universe with Einstein cosmology 
constant A. The metric of such local inertial systems is M-Beltrami metric that services as 
the basic metric of de Sitter invariant special relativity (dS-SR); In section 4, we derive the 
electric Coulomb law at light-cone of FRW Universe in terms of dS-SR Maxwell equations. 
As is well known, Coulomb force dominates the dynamics of the atomic spectrums. In section 
5, we discuss the fine-structure constant variation along the best-fit dipole direction shown 
in m- The a-varying Aa/ao = («z(D) — ao)/ao (where «o is cds value in laboratory) 
in this region is derived. Using the data along this best-fit dipole reported by [3111], the 
model’s parameters are determined. The theoretical predictions are consistent with the 
observations. In section 6, we examine the a-varying in whole sky. The results are also in 
agreement with the estimate from Keck- and VLT data. Finally, we briefly summarize and 
discuss our results. In Appendix A we briefly recall the Betrami metric and the de Sitter 
invariant special relativity. In Appendix B, a remark on the modified Beltrami metric used 
in this paper is provided. 

2 Light-Cone of Friedmann-Robertson-Walker Universe 

The isotropic and homogeneous cosmology solution of Einstein equation in GR (General 
Relativity) is Friedmann-Robertson-Walker (FRW) metric. In this section we discuss the 
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Light-Cone of FRW Universe because all visible quasars in sky must be located on it (see 
Figure [T]). 



Figure 1: Sketch of the light cone of the Friedmann-Robertson-Walker Universe. Only 3 coordinate axes 
{Q° = cf, Q 1 , Q 2 } are shown in this three dimensional figure. The Q 3 axis could be imagined. The Earth 
is located in the origin. The position vector for nucleus of atom between the QSO and the Earth is Q, and 
for electron is L. The distance between nucleus and electron is f ~ |L — Q |- The location of the Minkowski 
point of Betrami metric is denoted by notation “ x ” with M = (M°, M 1 , 0, 0). 


The Friedmann-Robertson-Walker (FRW) metric is (see, e.g., 

dr 2 


ds = cdt — a{ty 


1 — kr 2 


+ r dd + r sin 9d<j) 


= (dQ 0 ) 2 - a(t) 2 {dQ l dQ l + 


i \2 


k(Q l dQ 
1 - kQ i Q i 

= gAQ) d Q^QY ( 2 . 15 ) 

where a(t) is scale (or expansion) factor and r = yjQdQ 1 = Q, Q 1 = Q sin 9 cos 0, Q 2 = 
Q sin 9 sin Q 3 = Q cos 6 and (Q 0 ) 2 = c 2 t 2 . Hereafter, for the sake of convenience, we take 
t to be looking-back cosmologic time, so that t < 0. As is well know FRW metric satisfies 
homogeneity and isotropy principle of present day cosmology. For simplicity, we take k — 0 
and a{t) = 1/(1 + z(t)) (i.e., a(t o) = 1). And the red shift function z is determined by 


ACDM model 


where 


see, e.g., Eq.(64) of [29]): 

, , f° dz' 
t(z) = 


H(z')(l + z'Y 

H(z') = Ho\/^mo(^ + ^) 3 + ^2ro(1 + Z ') A + 1 
H 0 = 100 h ~ 100 x 0.705 km • s -1 /Mpc, 

~ 0.274, Q ro ~ 10" 5 . 


(2.16) 


n 


mOj 


a 


mO 


(2.17) 


Figure of t(z) of Eq. fl2.16l) is shown in Figure [2] 

The Light-Cone of FRW Universe is defined by ds 2 = 0. From Eq. fl2.15H . we have the 
light-cone equation: 

(dQ 0 ) 2 — a(t) 2 (dQ) 2 = 0, or — cdt — a(t)dQ — - — ^——dQ. (2-18) 

1 z(t) 
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-t(Gyr) 



Figure 2: The t — z relation in ACDM model lea. (12.161) 1. 


Substituting (12.161) into (12.18[) gives 

= (219) 
Figure of Q(z) of Eq. (l2.19|) is shown in figure [31 Ratio of Q over Q° is shown in figured 


Q(z) (Glyr) 



Figure 3: Function Q(z) in ACDM model lea. (12. 191) 1. 


o/(-e°) 



Figure 4: Function of Q(z)/Q°(z). Q(z ) and Q°(z) = ct are given in Eqs. (12.191) and (12.161) . 


3 Local Inertial Coordinate System in Light-Cone of 
FRW Universe with Einstein Cosmology Constant 

In principle, almost all calculations on quantum spectrums in atomic physics are achieved 
in the inertial coordinate systems. From the cosmological point of view, the phenomena 
of atomic spectrums should be described in the local inertial coordinate systems of FRW 
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Universe. Therefore, we are interested in how to determine the local inertial coordinate 
system in light-cone of FRW Universe when the Einstein cosmology constant A is present. 

Existence of local inertial coordinate system is required by the Equivalence Principle. The 
principle states that experiments in a sufficiently small falling laboratory, over a sufficiently 
short time, give results that are indistinguishable from those of the same experiments in 
an inertial frame in empty space of special relativity m Such a sufficiently small falling 
laboratory, over a sufficiently short time represents a local inertial coordinates system. This 
principle suggests that the local properties of curved spacetime should be indistinguishable 
from those of the spacetime with inertial metric of special relativity. A concrete expression 
of this ideal is the requirement that, given a metric g a y in one system of coordinates x a , at 
each point P of spacetime it is possible to introduce new coordinates x' a such that 


9ap( x 'p) = inertial metric of SR at x P , 

and the connection at x' P is the Christoffel symbols deduced from g' a p(x' P ). 

In usual Einstein’s general relativity (without A), the above expression is 


(3.20) 


'l',;,3< x 'ri = 1 h/3, and iV, = 0, 


(3.21) 


which satisfies the Einstein equation of E-GR in empty space: G lw = 0. 

In dS-GR (GR with a A), the local inertial coordinate system at x' P is characterized by 


9' a p{ x 'p) ~ Bap\ x 'p) — 




(M)( x / p 


+ 


a 


Vp\( x ~ M x )r] up (x P - M p ) 
R 2 cr( M )(x' P ) 2 


(3.22) 


with ^< )(x . p) = l _<hML-M^ P 


M v 


R 2 


a/3 


= + dpBW - aX"’) 


(M)1 


7T (8*r] up + S*r]p P )(x ,p - M p ), 


(or 


R 2 a( M \x' p ) M 
dg'^'p) _ x'" P - M v 


dx' P 


R 2 a( M \x' P ) 2 


‘^ j gpfi'9ay U gpoRBia T 9fio9pu T 


(3.23) 

^Vp^ypv'naxixp - M p )(x'£ - M x y 
R 2 <j( M \x' P ) 


where B^\x' P ) were given in fll.lOp . which satisfies the Einstein equation of dS-GR in empty 
spacetime: G pu + A g pu = 0 with A = 3/R 2 . (Note rj /iu does not satisfy that equation, i.e., 
+ A a/a 0. So it cannot be the metric of the local inertial system in dS-GR with 
A). 

To the light cone of FRW Universe with A, the coordinate-components Q°(z ) = ct(z ) 
and Q(z) have been shown in Eqs. fl2.16p (or Figure [2j) and (12.191) (or Figure [3]) respec¬ 
tively. Therefore from (I1.10|) . the space-time metric of the local inertial coordinate system 
at position Q(z) of the light cone is determined to be 


where 


B^\Q) = Bp V {Q-M) 


Vpu Vp\(Q X ~ M x )g up (Q p - M p ) 
cr( M )(Q) + R 2 a( M /Q) 2 


VpviQ 1 * ~ M P )(Q U - M v ) 
R 2 


(3.24) 


u (M) (Q) = a(Q-M) = 1 


(3.25) 






















We see from Figure |T] that the visible atom is embedded into the light cone at Q-point. Since 
Q ~ L (i.e., comparing with the Universe, atoms are very very small), we can reasonably 
treat the metric of the spacetime in the atomic region as a constant. This is just the adiabatic 
approximation adopted in [2l] . When Q M —>■ M M , we have B ( fl jP(Q) => r] pu . So, Q p = is 
the Minkowski point of the Beltromi metric b\^J\Q). 

Now let’s derive e p and u al ) from fid.241) . Setting 

q p = Q p - (3.26) 

then equations (13.241) and (13.2511 become 

R (M) _ q X Vu P q p , 97 x 

^ a (M) RZ^M)) 2 ’ [ J 

where 

ff (M) = ff (M )( ) = J _ (3.28) 

R z 

We introduce notations: 


% = q x = q x = 

t = v^q" = f M", 


and construct two project operators in spacetime {q 11 } with metric q pv \ 




q t JB 


Wv ,, _ 

—9 1 — _ 9 

qz gz 


It is easy to check the calculation rules for project operators: 

9 p x0 X v = 9 p x'q Xp 9 piJ = 9 pu , or in short 6-6 = 9, 
uj^xcD* v = oJ p \r] Xp co pu = uj pu , or in short cu • u = u, 
6 p xUJ X v = 6 p x6 Xp 6J P u = 0, or in short 9 ■ lu = 0, 

0 piJ + Cu pu = q pu , or in short 9 + u> = I. 

B { ™ ] (Q) can be written as 

g(M) _ _|_ q^qv 


a (M) R 2( a (M)}2> 


where 


cr 


(M) _ 1 


r\ i zvq p q v _ 1 _ _ 1 _ <?_ 


R 2 


R 2 


R 2 


(3.29) 

(3.30) 


(3.31) 


(3.32) 

(3.33) 

(3.34) 

(3.35) 


(3.36) 


(3.37) 

Since B p X J } is a tensor in the spacetime {g M , it can be written as follows from (I3.36jl : 

(3.38) 

), the above 


r?(M) _ \ n _i_ ^ 


Furthermore, by means of {B (My ) lllJ B^ X p = 8 P = q p and the rules (13.3111 - 
expression of B^ leads to: 

( j B {M)yv _ a {M)Q^ + 


(3.39) 
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Or explicitly in matrix form: 


/aW( 1-iff 


{( S (M))^} 


- g V<T (M) 

R 2 

-q 2 q 0 a( M '> 

W- 

- q 3 q°a <- M ) 
R 5 


-cr 


R 2 

<"»(! + 

_ ? 2 ? i CT (M) 

R 5 

- ? 3 gV M) 

R2 


In the Beltrami spacetime 15 with the metric B^\ the tetrad e“ A( is dehned via the following 
equation 

= W\e\. (3.41) 

Generally, we have expansion of e a fl : 

e\ = a6\ + boo\, (3.42) 

where a and b are unknown constants. Substituting equations (I3.38P (I3.42p into (13.411) gives 


_q 0q 2 a (M) 

R 2 


^ql q 2 a (M) 


R? 


_ q 0 q 3 a (M) 

R 2 

_ g l g 3 ff (M) 


\ 


-^(1 + T 

_ q 3q2 (T (M) 

R 2 


( 9 2 ) 2 


R 2 

_ q 2 q 3 a (M) 

i + V )/ 


(3.40) 


®uv + 


a 


(AT) p v 


(cr 


(M)) 


2 <*V = W a 0° tt + 6w °u)( a 0 6 i/ + 6t *O 


= a e,„, + 6 


/XI/* 


(3.43) 


Comparing the left side of (13.431) with the right side, and noting 6 and a; being project 
operators with properties of (I3.32D - 03.35I) . we find that: 


a = 


a 


{MY 


b = 


a 


(M)' 


(3.44) 


Substituting (I3.44D into (I3.42p gives 


e R = 


a 


;0 a „ + 


co 


(M) /4 1 fj(M) M 


1 \ ^(g A -M A )(Q"-M") 


cr(M)> 1 \ a (M) 

e 0 M is the inverse of e“ given by: 


(M) 


(1 - a( M ))i? 2 


(3.45) 


e ' , e° = S a 

C 'a ° (i u a ’ 


(3.46) 


With the equations (I3.46P and (I3.4ip . we have 

= + 

, ° (M) - ri ab s b x(Q x - m*)(q» - m») 
a+ 1-aM R 2 


(3.47) 
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and 


_ 


= e\(B< M >r = {\J^ m S a „ + 2jy<S“ Ma <M> S"» + (a< M >Yu,-“) 


= Va^9 afl + a {M) Cu ai1 

= ../y^y i <rlM> ~ ^ ^ - m " ) 2 (q " ~ M “\ 


Cdfi 


9 an + 


a (M) 1 (j(M) 


OJ, 


afi 


( 1 1 V«a^(Q A -M A )(Q^-M") 

V a ( M ) a?7i//i V ff(M) (1 — a( M ))i ? 2 

Next we derive spin-connection cu a? ) t . From identity 


e a ;j/ — ^ e a + v e b + — 0, 

1 

Xfl ~ 2 


rL = kB {M) T(dxBM + S„B<f > - 


we have 


<. = - e^) - ir$„(e“VJ - e“ej). 


'P' 2 ^A tV “ “P 

Substituting (j3.27[) (I3.45P into (I3.5ip and (I3.52p gives 

1 


oj ab = 
p 


R 2 ^1 + \/VaW 


(«-«)(g A -M A ). 


(3.48) 

(3.49) 

(3.50) 

(3.51) 

(3.52) 

(3.53) 


4 Electric Coulomb Law at Light-Cone of FRW Uni¬ 
verse 

The hydrogen atom is a bound state of proton and electron. The electric Coulomb potential 
binds them together. The action for deriving that potential of proton located at Q = Q M = 
{Q° = ct , Q 1 , Q 2 , Q 3 } with background space-time metric g liu = b\!^\Q) of eq. f|3.24p ( see 
FigiB in the Gaussian system of units reads 

5 = / d 4 LV^F IJV F>- - e - J d^L^gfA,, (4.54) 

where g = det (B^), F, xv = and j' Jl = {j° = cp proton / \JB { 0 J /) , j} is the 4-current 

density vector of proton (see, e.g, Ref. [38]: Chapter 4; Chapter 10, Eq.(90.3)). Making 
space-time variable change of // — > LA — = {x° = — ct, x l = L l — Q 1 } and 

noting ~ we have action S' as 

J d‘xC^(B { F(Q))F„ u F^ - e - j i‘x\J-&et(B { ™\Q))fA^ 

= (OB^’tC) J -) / J-Aet(B$P(Q )), 

(4.55) 
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and the equation of motion SS/SA^x) = 0 as follows (see, e.g, [38], Eq. (90.6), pp25T) 


Ajr 

d v F> M ' = (B^ M) Y x d„F\ = ~—f- 


(4.56) 


In Beltrami space, A M = {4>b, A} (see, e.g., [38], eq.(16.2) in pp. ^5) and 4-charge current 

= {cpproton/ \JB^\ j}. According to the expression of charge density in curved space in 
Ref. [38], ( pp.256 , Eq. (904)), /Voton = Pb = ^=<5 (3) ( x ) and j = 0, where 

7 = det( 7 ij), (4.57) 


dl 2 = 7 ijdx'dx? = —B): + 


r>(M) td{M)' 

(M) , n 0i n j 0 


B. 


(M) 


dx'dx 1 (see eg.(84.7) in Ref.[ 19]04.58) 


oo 


1. When p — i (7 = 1,2,3) in Eq. ('14.56]) , we have 

A r-rr 

d l d tl A 11 ~ / = 0 . 

By means of the gauge condition 


(4.59) 

(4.60) 

(4.61) 

A* = 0 (4.62) 

is a solution that satisfies the gauge condition (15.79[) (noting d 0 A° = 0 b (tb) = 0 

due to *£ - - 


we have 

Then 


= 0 , 

(B^y^d^A* = o. 


-ggy — -Tfij = 0 ). Eq. (14.621) is the vector potential. 

2. When p = 0 in Eq. fj4.56p . we have the Coulomb’s law: 

(B^f(Q)d i d,Mx) = = 

C C / r->(M) 


- 47 re 


B, 


—4ne 


oo 

<S (3) (x), 




:<S <3| 0 


(4.63) 


det (B$\Q)) 

where B^'y = —det(5^) has been used, and b\^ ] (and ( B ( M ))b) were gi ve n in 
(1^271) (and (137401) ). he., 


{(bW)«) = 


( -^( 1 + ^ft 


V 


-gV<7 M > 

R 2 

-, 3 oh (M ) 

R 2 


—q 1 q 2 g( M ) 

f+¥) 

_q3 q 2 a (M) 

R 2 


a^ M \ 1 + ^f) / 


_ ? i g 3 0 .(M) 

R 2 

_ q 2 q 3 a (M) 

R 2 
(9 3 ) 2 




(4.64) 


The equation of Coulomb law (j4.63j) can be compactly rewritten as: 

- (V^ (M ) Va; ) = -rA 

where tensor 93 := {(£h M ))b\, operator V := {A}, 77 = . == and 

\/det(A„ J (Q)) 


(4.65) 
3 := 

<A)(x). Symmetric matrix 93 (iU ) can be diagonalized via similitude transformation 
due to matrix P: 

Ai 0 0 


p T^(M) p = Ab= [ 0 A 2 0 

0 0 A 3 


(4.66) 


12 





































Here, A* with i — 1, 2, 3 and matrix P can be found in det(23^ M ) — AI) = 0. From 
(I4.64p . the results are: 

(q) 2 + R 2 


A, = -a^- 


R 2 


A 2 = A 3 = -cr (M) , 


P = 


( P\/(<1 3 ) 2 _ q 1 q 2 _ g 3 

q 3 yRf_ A / q 2 [(g 1 ) 2 + (g 3 ) 2 ] \/ (<7 1 ) 2 + (q 3 ) 2 

g 2 yW \/ (g 1 ) 2 +(g 3 ) 2 




V 


yj (9 3 ) 2 


jV 


0 

n/IT ) 2 


(4.67) 


(4.68) 


Vq 2 [(9 1 ) 2 + (9 3 ) 2 ] V(9 1 ) 2 +(9 3 ) 2 

where q* = Q l — M l with % = {1, 2, 3} (see Eq. (l3.26jl ). and q 2 = (q 1 ) 2 + (q 2 ) 2 + (q 3 ) 2 . 
The follows can be checked: 

P T P = PP T = I, (4.69) 

where / is 3 x 3-unit matrix. So that Eg. (14. 65 j) can be rewritten as follows 

- (V T X PP T & M) PP T V X ) (j)B — (( y T x P)(P T ^ M) P)(P T V x )) fa 

= - (VJ’AbVJ 0b = -y3 = 7 ~ 47re ,„ =^ ( 3 ) (x), (4.70) 


det( 4 f ) (Q)) 


where 


y = Px, or y l = P^, x l = P^y ] 


V, y = P T V X , or d/dy 1 = Pj,dfdx\ d/dx l = Pijd/dy j . 


(4.71) 

(4.72) 


Substituting Eqs. fl4.66p (I4.7ip into Eq. (l4.70p gives 

d 2 a 2 a 2 

+ WhrW + (v/v^v) 2 . 


<Pb = -r lS(l%v i )6(F%vi)6(F%y i ) 


= -V 


S{y l )5{y 2 W 

|det(if)| 


- 47 re 




A 2 2 


-A 3 ' 


det(4f } (Q)) 


\/— A 1 A 2 A 3 


(4.73) 


where a calculation result |det(P^)| = 1 has been used. Setting 


x = 


y 




= with i — 1, 2, 3, 


(4.74) 


(I4.73p becomes 


d 2 a 2 
+ 


a 2 


= vit 


^(x 1 ) 2 d(x 2 ) 2 <9(x 3 ) 2 

—47re <5(x 1 )5(x 2 )(5(x 3 ) 


det(S^ ) (Q)) 

Then, we get the Coulomb potential: 

1 


a/—A 1 A 2 A 3 


(4.75) 


0 B — 0 b(Q) — 


yJ-det(B^\Q))^~ XlX 2 Xsr B : 

where = ^/(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 . 


(4.76) 
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5 Fine Structure Constant Variation in {Q°, Q\ o, 0 } 

At this stage, in order to get the expression of fine-structure constant a(Q) = ct 2 (f2) 
at point Q = {Q° : Q 1 , Q 2 1 Q 3 }, we should substitute Equations (I3.47|h (I3.53|) . (I4.62p . 
(I4.76P and (11.10p for e a M , oj ab , A\ A 0 = cj)B and Bj^' respectively, into the Dirac equation 
(1 1.12 p of hydrogen atom located in the local inertial coordinate system of the light-cone in 
FRW Universe with A (see Fig{l]). Such an a(Q) will characterize the temporal and spatial 
variations of fine-structure constant a. When we assume M — 0, a(Q)|(M=o) has been 
calculated in the Ref. [211122] • In the present paper we study the observation results of Keck 
and VLT reported by |3j recently by means of M = {M°, M 1 , 0, 0}-model. 

To calculate a(Q) analytically, we build 3-dimension spatial Cartesian coordinate sys¬ 
tem {Q 1 , Q 2 , Q 3 } on the equatorial coordinate figure showing the Keck and VLT best-fit 
dipole structure of A a/a in Ref. [3j. Denoting Q l = Q*/|Q| and Q 1 as the best-fit dipole 
position, the directions of the three axis {Q 1 , Q 2 , Q 3 } hi this figure are t?^) = 

{17.4 h, —59°}, Q 2 ((pf A \d Jp) = {17.4 h, 31°}, Q 3 {^ A \ = {11.4 h, 0°} (see Figure 

[5] and its caption). Angle (0) between a quasar sight line ({cp]/ 64 ', and axis Q 1 is 

determined by following 

[RA] _ [RA] 

cos 0 = cospdf^] cos[$J^] cos[—-- — tt\ + sin[i?^] sinft?^]. (5.77) 

In this section we calculate ct(Q) for 0 = tt and 0 = 0 following the method in Ref. pHlj 
and ||22]. 



Figure 5: (color online) The 3-dimension spatial Cartesian coordinate frame {Q 1 , Q 2 , Q 3 } on the equatorial 
coordinates. In left figure (a), the background is all-sky plot in equatorial coordinates showing the indepen¬ 
dent Keck (green, leftmost) and VLT (blue, rightmost) best-fit dipoles, and the combined samples (red, cen¬ 
ter), for the dipole model Aa/a = dcos© copied from [3j. The locations of axis {Q 1 , Q 2 , Q 3 } in this figure 
are marked by . In [3] it has been measured that the best-fit dipole is at right ascension tp{ RA 1 = 17.4±0.9 
h, declination = —58 ± 9 dec. The cosmic microwave background dipole antipole are illustrated for 
comparison. The directions of 3-dimension spatial Cartesian coordinate system {Q 1 , Q 2 1 Q 3 } that we take 


are: Q 1 (y 


[RA] ADU _ 


1 > ^1 


) = {17.4 h, -59°}, Q 2 (y 2 


2 f, n [RA] ADU _ 


< J ) = {17.4h, 31°}, 0 3 (y : 


3/.ARA] ADU _ 


3 )$3 


) = {11.4 h, 0°}- In 


right figure ( b ), the 3-dimension spatial Cartesian coordinate system {Q 1 , Q 2 , Q 3 } with a non-zero space 
component M 1 of Minkowski point is drawn. 0 is angle between a quasar sight line and axis 

Q 1 . Formula for computing it is cos© = cos[d] D ^] cospd^ 1 


cos 
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5.1 Formulation of Alp ha-Variation for Case of 0 = 7r, and 0 

When quasar sight line is anti-parallel (or parallel) to direction of Q x {tp'\ iA \ d, 7 ^) = 
{17.4 h, —59°}, the angle 0 between them is equal to n (or 0), and the locations of distant 
atoms on the light-cone in the FRW Universe are at {Q°, Q 1 > 0 (or < 0), Q 2 = 0, Q 3 = 0} 
(see Figured]). Then we have q° = Q° — M°, q 1 = Q 1 — M 1 , q 2 = 0, q 3 = 0, and 


b$HQ) 


/ i , (Q°—M 0 ) 2 
I" Jf(,(M])2 
(q°-M°)(Q 1 -M 1 ) 
R 2 {ct(M)) 2 

0 

v° 


(Q°-M°)(Q 1 -M 1 ) 
R 2 {a( M )) 2 
-1 , {Q 1 -M 1 ) 2 

-JMJ -I" #2( ct (M))2 

0 

0 


0 

0 


Cj(M) 

0 


0 ^ 
0 
0 

/ 


(5.78) 


( g(M)y"(Q) 


with 


/ 

V 


(Q°—M°)(Q 1 —M 1 )ct ( m) 
R 2 


0 

0 


(Q°-M°)(Q 1 -M 1 )c r< M ) 
R 2 


0 

0 


0 0 \ 

0 0 

0 

0 -aW / 


(5.79) 


a (M) 


(Q° - M 0 ) 2 - (g 1 - M 1 ) 2 

7u 


(5.80) 


Here, for a known red-shift z, g° = ct(z) and Q l = ±^/g(z ) 2 — (Q 2 ) 2 — (Q 3 ) 2 = ±Q(z) 
are determined from Eqs. fl2.16p and (12.19[) respectively (see also Fig ]2] and Figj3]). Then 
Eqs. fl4.67h -f l4.68h become 


A, = -V«> (l + 1)2 ) , A 2 = A 3 = (5.81) 

1 0 0 \ 

0 10, (5.82) 

0 0 1 / 



Substituting f!3.47p f!3.53jl fl4.62l) (14.76(1 and fjl. 10(1 into (11. 12(1 gives dS-SR Dirac equation for 
the electron in the distant Hydrogen located at the light-cone of FRW Universe: 


hc/3 


._ n (M) _ J~Z{M) rn r 

V L U + i — - VabS a x(Q X - M A ) 7 b (g^ - 


(1 - a^)R 2 


h 


fp = o, 


(5.83) 

where factor hc/3 in the front of the equation is only for convenience, U 1 ~ Q M has been used 
(see Figjl]), and 

v L = -3- 

d i 1 

-'V 


8L, 


dL “ 4 i? 2 (l + VJW 


(i‘A - - M x )o at 


,. e 5„„B'" l (S)+Vi-Sffl” , (S) 


(«), 


ch 


det (B$P(Q)) 


y/— AiA2A3 Tb 


(5.84) 


We use the method suggested in EU 122 ! and expand each terms of (15.831) to the order as 
follows: 
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1. Since observed distance hydrogen atom must be in the light cone and the location is 
{Q°, Q 1 , Q 2 = 0, Q 3 = 0}, then rj ab L a L b ~ r] ab Q a Q b = ( Q °) 2 - (Q 1 ) 2 , and the first 
term of (15.83j) reads 


hcpiV^Yvj; = \Ji- ^ Q0 M °) 2 r2 ^ Q1 (5. 


with PY = {^ 7 ° = p 2 = 1, py = a*} 


85) 


(5.86) 


c) fie 3 

hcfiiyvfa = (ift— + ihea ■ V + —^-ufya ab + ePYB llv A v )ip, (5.87) 
In the follows, we use variables {x 1 , x 2 , x 3 } (where since Eqs. (l5.8ip(l5.82p . x l = x l / v 7 — X\ 


we have x 1 = 


1 ~2 _ 1 


[1+(Q 1 -M 1 ) 2 /R 2 ) 


x , x = 


y/ <j( M ) 


x 2 and x 3 = 


3 — _i 




{#,£ 2 , x 3 }. Following notations are introduced hereafter: 

r B = ix 1 + jx 2 + kx 3 , \r B \=r B , 

d . d , <9 

J— + k 


Vr — i 


dx 1 dx 2 


clx 3 ’ 


x* G {x 1 , x 2 , x 3 }. 


x 3 ) to replace 

(5.88) 

(5.89) 


d d 


Then, noting Qxl Qxl Q£l ^ ct( m)( 1 + (qi_mi) 2 /.R 2 ) dx 1 > ax' 
i = {2, 3}, the eq. 05.87p becomes 

hcfti^'D^ip = ( ih-fr —t i—?==d ■ Vg + i 


<9x z d 
dx x dx 1 


1 a 






vTPO dxi 

i d 


with 


s/l + {Q 1 - 


- 1 


a 


dx 1 


. ^ C /3 ab y. , d (M) 

H 7 ^ Cab + eB 00 


4>b + ea 1 b[q I] < f> B ^j 


( d he _ he 

= in— -f -1 , a ■ \7 b + » . 

I dtL \J cr( M ) v 7 (j( M ) 


■\/l + (Q 1 - M l ) 2 /R 2 


- 1 


a 


di 1 


+ ^<7^ab + 


1 


(Q° - M°) 

ct( m ) 1 i? 2 (cr( M )) 2 


21 


+ 


e0s — 


(Q 1 - M 1 )(Q° - M°) j 


l ? 2 (( j ( M )) 2 


a ecj>B ) ' 0 - 


(5.90) 


2. Estimation of the contributions of the fourth term in RSH of (I5.90p (the spin-connection 
contributions) is as follows: From (13 . 531) . the ratio of the fourth term to the first term 
of (157901) is: 

hy^a.b^Ufjab^ he Ct 1 Ct k 1 ct(l r 


ihdtip 


4 2 R 2 rrifC 2 8R 2 m e c 8 R 2 


0 , 


(5.91) 


where a c = h/(m e c ) ~ 0.3xl0” 12 m is the Compton wave length of electron. O(cta 0 /R 2 )- 
term is neglectable. Therefore the 3-rd term in RSH of (15.871) has no contribution to 
our approximation calculations. 


3. Substituting (I5.9ip (I5.90P (I5.86P into (I5.85P and noting L a ~ Q a (see FIG. 1), we get 
the first term in LHS of (I5.83p 

he 


/ Q 

hc/3iVcrWyvfy = VffW 1 ih— -b i 

V o£.L 


+i 


he 


1 


yj\ + (Q 1 -M 1 ) 2 /R 2 


- 1 


v 7 <j( M ) 
i ^ 


a ■ Vr 


a 


dx 1 


(5,92) 


_a( M ) R 2 (eH M )) 


R 2 (a( M y 
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4. The second term of (15.831) is 


— -\J ff(M) 

hcpi- - Vab(Q a - M“) 7 b (Q^ - M»)Vfy 


(l-o-( M ))i? 2 

cr( M ) - r 


= Hcfii 


7 °(Q° - M°) - j ■ (Q - M) [(Q° - M°)V 0 + (Q i - M l )D t \ p 


(1 — cr( M ))i?2 L' 
a {M) - vV(M) 

[(Q ~ M ) ~ & {Q ~ M )\ 

7°7 1 (Q 1 — Tf 1 ) _ i _ i _ 

2it 2 (l + V^W)v^(My “ + i? 2( (T (M))2 

7 ° 7 1 (Q°-M 0 ) *e (Q° - M°)(Q 1 - M 1 ) 


(Q° - M°) ^ - 

+(Q 1 -M 1 ) (V + 


*e 1 


(Q° — M 0 ) 2 


2RR1 + ch i? 2 (cr( M )) 2 


V' 


ct (M) _ 

(1 — a( M ))i ? 2 
x < iTic 


[((/-M^-cd^-M 1 )] 


<«° - M0) aF + w ‘- Ml) aF 


g°-M» ^ , , 
( ct (m ))2 e< ^ B r"^> 


(5.93) 


where Q 2 = Q 3 = M 2 = M 3 = 0, 7 0 = /3, /Ty 1 = a 1 and Eqs. (l5.84j) (I5.80j) were used. 

Nntino- 'rM = />„/> d = -4L = 9 9 — _JL _ 9x)__9_ _ 1 9 

° ’ 9L° &e° c9ti, ’ 5L 1 cte 1 Si 1 9a; 1 9a; 1 > 

(I5.93|) becomes 

fiqfl (1 _ (j(M))j^2 W5(Q A - ^ VW - M*)T%1> 

crW _ v/TiMf 


X 


(g° - + (g 1 - m 1 ) ^ 

dt L W v/a( M )(l + (g 1 - M 1 ) 2 /R 2 ) dx 1 

, Q°-M° , ) , 

+ ( a (M)) 2 e ^J^ 


(5.94) 


5. Therefore, substituting (15.92p (I5.94p into (15.831) . we have 


, ,_ a {M) _ J a (M) 

ih ( + —-^ 7 TTTvr[(g° - M°f - (g° - M°)(g 1 - M^a 1 } 


= < —ihca ■ V b 


—ihc 


(1 - a( M ))R 2 

1 


d 

dtr 




+ 


(g° - M 0 ) 2 (j( M ) - (go - M 0 ) 2 


+ 


Va( M ) R 2 a^ M Wa^ M ) (1 - <x (M) )R 2 (cd M )) 2 

1 


ecj) B + m e c 2 f3 > i> 




, , i 9 , (g 1 -M 1 )(g°-M°) , ^ 

— 1 | a r H-:——— , - a ecps 


dx 1 


a m _ 

(i - V^))r 2 


ihc{Q 1 Q° - (g 1 - M 1 ) 2 ^ 1 ) <9 
vV M )(l + (g 1 - M 1 ) 2 / R 2 )'d5d 


R 2 a {M)yfcZM) 


fp- 


(5.95) 


In order to discuss the spectra of Hydrogen atom in the dS-SR Dirac equation, we need 
to End its solutions with certain energy E for electron in the atom. From Eq.( 68 ) in 
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Ref. m • we have 


p° = — = ih 
c 


(Q° - M°)(Q U - M u )\ a 
1 R 2 ' " 


5(Q° - M°) 

2 R 2 


E = ih 


(Q° — M 0 ) 2 c(Q° — M°)(Q 1 — M 1 ) 5 c(Q°-M°) 

-^-7V7-+ 


R 2 


2R 2 


Eip~ih[i- w -m, cW ° - M ° )W ~ Ml) 


R 2 


d 


R 2 + (Q 1 - M 1 ) 2 //? 2 dx 1 


(5.96) 


where a estimation for the ratio of the 4-th term to the 2-nd of E'lp were used: 


\ih^\ 

I — c 2 t 2 


I Ul 2R? I 


5 h 


5 a. 


~W 


ihd tL i>\ 


-2 c 2 t 2 
R 2 


E I 


2 tmpC 2 2 ct 


where a c = h/(m e c ) ~ 0.3 x 10 _12 m is the Compton wave length of electron and ct is 
about the distance between earth and a distant atom near quasar. Obviously, a c /(ct ) 
is ignorable. For instance, to a atom with ct ~ 10 9 ly, a c /(ct ) ~ 10 ~ 38 << (ct) 2 / R 2 ~ 
10” 5 . Hence the 4-th term of E'lfj were ignored. Eq. fl5.96p means 


d 


= T 
dt L 1 


E 




- (Q° - M°) 2 /R 2 
ihc(Q l - M 1 )(Q° - M°) 


d 


R 2 (l ~ (Q° - M 0 ) 2 /R 2 ) v /<j ( M )( 1 + (Q 1 — M l ) 2 /R 2 ) dx 1 

Then substituting (I5.97[) into (I5.95[) gives 

Et/j = Hoip + H'xj) 


^•(5.97) 


(5.98) 


with Hn =11 — 


(Q° - M 0 ) 2 


R 2 

x —ihea ■ Vb - 


vVW 


- Vcr^)(Q 0 - M°f 


(l-cr( M ))i? 2 

(Q° — M 0 ) 2 ( 1 


-1 


i?V M ) + 


ecf>B + m e c 2 p 


e ffO 2 

= —ih z (£L)ca - Vs--- 1 - m e , z(tl)c 2 /3, 


tb 


(5.99) 


H'= 1 - 


(Q° - M°) ; 
tf 2 


1 


( ct (M) _ Vcr( M ))(g° - M 0 ) 2 


-1 


v /1 + (g 1 - M 1 ) 2 /# 2 1 J “ dx 1 


(*((g x - M 1 )(g° - m°) - (g 1 - M 1 ) 2 a 1 )^ T - 


x ^ — */ic 

a (M) _ 

' (1 -cr( M ))R 2 

(cd M ) - v / ^))(g° - M 0 )^ 1 - M 1 )a 1 

+ (i - (g° - m°) 2 /^ 2 )(i - cr (M) )^ 2 

ih^Q 1 - M 1 )(Q° - M°) d 
~ i?Ved M Hl + (Q 1 ^ ERY/R 2 ) di 1 

3 

= 


(1 _ a ( M ))R2 

(g 1 -M 1 )(g°-M°) i 


__ a ed>B 

R 2 (J (M) v ^(mT 

i. d (g 1 -M 1 )(g°-M 0 ) 2 


1 £ + i/ic 


(<j ( M )) 2 

(Q° - M°)(Q 1 - M 1 ) d' 
R 2 y/1 + (g 1 - M 1 ) 2 /R 2 di 1 




(5.100) 
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where 


and so that 


Q l = a\ C^Oil/R 2 )- Q 2 = -^ T , C 2 ocO( 1/R 2 )- 
Q 3 = a 1 -^, C 3 (X 0(1/R*), 


H' oc 0(1/R 2 ) « H 0 oc 0(1). 


(5.101) 


(5.102) 


We have mentioned in the Introduction section that the operator-structure of H 0 of 
(I5.99P makes the corresponding eigen-equation Eo'i/j = Hoijj to be integrable. Hence 
Eq. (15.1021) means that H 0 (Eq. fl5.99p ) and H' (Eq. (j5.100ft ) can be legally treated as un- 
perturbation Hamiltonian and perturbation Hamiltonian respectively in QM-problem 
with H = H 0 + H'. 


6 . From Eq. fl5.99p . we have 


h z (Q)c — He i 1 — 


(Q° - M°f 

I? 2 


V^j + 


(^(m) _ v^M)(go -M 0 ) 2 
(1 - a^)R 2 


-1 

,(5.103) 


= e 2 1 


x 


(Q° - M°f 

Jr 

(Q° - M°) 


./yvn | (a (M) - V^))(Q° - M") 


y^M)+ R2(j{M) \ vl + V ^(M)yJ ^ 




(1 - (j( M ))i? 2 

1 


-1 


det (B^ ) (Q)) 


\f— A1A2A3 f 


where (^-expression (I4.76P were used, and hence 


(Q°-M°) 2 


ay(fi) = a 


sRRM) ^ R 2 a( M ) (i +v 7(M)i e 


-det(4f } (g))V-AiA 2 A 3 


, with a — —. 

he 


(5.104) 


(5.105) 


For case of 0 = n, (or 0), B\JJ\Q) and Ai, A 2 , A 3 have been given in Eqs. fl5.78p and 
f!5.8ip respectively, and hence we have 


det (b\™\q)) 


o 


(M) 


1 (Q 1 - M 1 ) 2 

W> J R 2 ( C r(W)2 


1 + (Q°-M°) 2 \ | (Q° - M°) 2 (Q 1 - M 1 ) 2 


cr 


(M) R 2 (a ( ~ M '>)‘‘ 


i? 4 (cr( M )) 4 




(5.106) 

(5.107) 


Substituting (15.106ft (15.107P into (15.105P gives 

1 L (Q°—M 0 ) 2 ( 1 \ 


ay(Q) = a 


1 (Q 1 -M 1 ) 2 \ ( 1 , ( Q°-M°) 2 \ , (QO-MOpCQi-M 1 ) 2 

R 2 ( ct ( m )) 2 i l + R 2 (a( M )) 2 J + J R4( cr (A/))4 


1 + 


(Qi-M 1 ) 2 

R 2 


(5.108) 
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When z = 0 (or Q° —* 0, and Q 1 —* 0), ct 2 (f2) should be normalized to be a 2 (f2 )| 2=0 = 
a 0 which is the a-value measured in the Earth laboratories. So from (15.108jl we have 


AM) 

T o 


+ 


(M 0 ) 2 


oiq — a- 




AM) 


(M 1 ) 2 


_L_ _|_ 

(M'l l 


(M 0 ) 2 


R2^M)y J l a (M) ^ R 2 (a[ M) ) 2 / K 4 (4 M) ) 4 


+ 


(M°) 2 (M 1 ) 2 


1 + 


(M 1 ) 2 

/? 2 


= aiVn 


(5.109) 


where 


(M 0 ) 2 


An = 


^ + jp.s“>v^r ViV'S" 


(M 1 ) 2 


^ R 2 ^r ; y 


ITT + 


T 77 ' « 2 (-r ; ) 


( M °) 2 


IMK 


(M) (M) I 


aA = a 


= 1 


(PI 0 ) 2 - (PI 1 ) 


(M°) 2 (M 1 ) 2 

H 4 (4 M) )4 


=(5.110) 


1 , W 1 ) 2 

1 + WF" 


\ Q°=Qi=0 Jl 2 

Therefore, from (15. 108jl and (15. 109j) . we obtain 
Aa _ ay(f2) — a 0 


(5.111) 


«o 


«o 


(Q°-M °) 2 


cr(M) T R 2 a (M)y/^(M) l l+y/^(M) 


Nr, 


(Q 1 —M 1 ) 2 


a(M) R 2 (a(M)) 2 I \ct(m) t jR 2( (T (m))2 


1 , (Q°—M 0 ) 2 \ . (Q°—M 0 ) 2 (Q 1 —M 1 ) 2 

I D2 / \2 I I 


i? 4( cr (M))4 


1 + 


(Q 1 —M 1 ) 2 
R 2 

(5.112) 

= ±Q(© (i-e., 


where Q 0 ©) = c£(z) and Q x (z) = ± a/Q ©) 2 — (Q 2 ) 2 ~ ( Q 3 ) 2 

V ' _ (Q2 = Q3 =0 ) 

IQ 1 ©)! = Q( 2 ;)) have been given in Eqs. fl2.16p and (12.191) respectively (see also Figj2] 
and Figj3]). 


5.2 Comparisons with Observations of Alpha-varying for Cases of 

© = 0, 7T 

Equation (j5. 1 12j) is the prediction of a-varying derived from dS-SR Dirac equation of 
distant hydrogen atom located in Q 1 -axis. When the location Q 1 > 0, the direction of the 
corresponding quasar sight line is opposite to the direction of Q^axis (see Figj5j). For this 
case the corresponding angle ( 0 ) between the quasar sight line and the Q 1 -axis is equal to 
7 r (i.e., 0 = 7T for this case). Oppositely, when the location Q 1 < 0, we have 0 = 0, i.e., the 
direction of the quasar sight line is same to the direction of <5 1 -axis. For both of cases, we 
can generally write Q 1 in Eq. fl5.112p as 

Q 1 = — IQ 1 1 cos 0 = — Q(z) cos 0. (5.113) 

Looking back cosmic time variable is Q° = ct < 0, and the coordinates of Minkowski point 
of Beltrami metric are M° < 0, M 1 > 0, and M 2 = M 3 = 0. Substituting (15. 1 12)1 into 
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(15.1121) gives 


Aa 

«o 


1 , (Q°(z) — M 0 ) 2 

( 1 'll 


yi + VcrW J\ 


/N 0 


T 


(-Q(z) cos Q-M 1 ) 2 V 1 , (Q°(z)-M°) 2 \ 

fl 2 ( ff ( M ))2 + i J 2 ( CT ( M ))2 / 


- 1 , 


(Q°(z)-M () ) 2 (-Q(z)cose-M 1 ) 2 Y-, , (-Q(z) cos 0—M 1 ) 2 \ 

R*(<j <* 0)4 + 7? 2 / 

(5.114) 


where 7V 0 is the same as (15.11011 and 


= 1 - 


(Q°(z) — M 0 ) 2 


(- cosQQ(z) — M 1 ) 2 
R 2 


(5.115) 


The equation (15. 114ft is the theoretic prediction of dS-SR, whose variables are z (red shift) 
and ©(= 0, or 7r), and the adjustable parameters are R, M°, M 1 . The discussions on it are 
follows: 


1. Since \R\ is the maximal length scale parameter in the theory (say |R| ~ 10 12 lyr [25]). 
we could deduce the Acc/cto’s Taylor-power series of l/R 2 from (15.1141) (for practical 
calculations, “Mathematica” is useful): 

— ~ -^-r \(M 1 ) 2 Q°(z)(-2M° + Q°(z)) + 2M 1 (M° - Q°(z)) 2 Q(z ) cos© 
a o 8 it 4 

+ (M° - Q°(z)) 2 (Q(z)) 2 cos 2 0] + 0(1/R 6 ), (5.116) 

where the leading term oc 0(1/R 4 ) is dominating in the expansion of Eq. (15.114p . 
Suppose the parameters and the variables are chosen such that 

|M°| » | M\ Q°(z ) ~ Q(z) - e, (5.117) 

we have 

— ~Tj 2 co S 0 M'(M»)V (6.118) 

a 0 8 /f 4 

where the lesser terms oc O((M 1 ) 2 M 0 e/R 4 ) and oc 0(e 2 / R 4 ) have been ignored, and © 
is only to be 0 or n (and noting cos(0) = 1, cos( 7 r) = —1). Thus, when 0 = 0, we have 
Aa/cto > 0 , and when 0 = 7 r, oppositely, we have Act/cto < 0 . This is interesting since 
Eq. (15. 1181) indicates that the scenario reported by [3] could be interpreted by (15.114|) 
with a particular parameter setting in proper region of variables (15.1171) in the theory. 

2. In this scenario we need to determine the parameters R, M°, M 1 by comparing the 
theoretical predictions with the observation data. Keck+VLT data have shown the 
relations between Act/ao and r = ct(z) along Q 1 -axis [3|. Let’s use Eqs. (l5.114p and 
(15.115P with © = 7 r and 0 to fit the figure 3 in [3] which is based on the combined Keck 
and VLT data and t(z) expression (12.161) . The best fitting gives 

R = 500 Glyr = 0.5 x 10 12 Lyr, 

M° ~ -100 Glyr = -1.0 x 10 11 Lyr, (5.119) 

M 1 ~ -22 Glyr = -2.2 x 10 10 Lyr, 


which are consistent with requirement of (15 .1171) . The fitted curve of (I5.114p with 
(15.115P is shown in Fig [6] 
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Figure 6: Determination of parameters R, M°, M 1 via fitting Keck+VLT’s a-varying data reported by [3]. 
Aa/cto vs A r(z)cos& with 0 = {0, and 7r} shows an apparent gradient in a along the best-fit dipole. 
The best-fit direction is at Q 1 (ip[ R/ ^ , h{^) = {17.4 h, —59°} (see figE]). The data reported by [3] are 
shown with error bars. A spatial gradient is statistically preferred over a monopole-only model at the 4.2cr 
level. A cosmology with parameters (Ho, Dm, Da) were given in (12.171) . The fitted model’s parameters are 
R ~ 500GLyr, M° ~ — lOOGLyr, M 1 ~ —22GLyr. The resulting curve of Aa/ao(r(z)) of Eg. (15.11411 is 
shown. 


Table 1: Examples of predictions of Aa/oto'. 0 = {0, 7r} is angle between quasar sight line and axis Q 1 . 


For each redshift z, there is a pair of (Aa/cno)th-predictions from Ea. (15.1141) with parameters in Eo. (15. 1191). 


z 

0.65 

1.47 

2.84 

0 

0 

71 

0 

71 

0 

71 

(Aa/a 0 )th 

0.60 x 10” 5 

-0.62 x 10~ 5 

1.20 x 10" 5 

-0.88 x 10" 5 

1.88 x 10~ 5 

-0.98 x 10”6 


3. After determination of R , M°, M 1 of Eq. (l5.119p . the metric B^J\Q) of local inertial 
coordinate system and then Aa/ao(z) along Q 1 -axis are fully known. In figure [3 
the curves of A a/ao(z) are plotted. The Keck’s data in 2004 reported by [9] [39] are 
illustrated for comparison. Since the 2004-data were obtained by observations in all 
directions in Keck at that time, the deviations between the data and the prediction 
curves of Aa/ao(z ) are understandable. The point here is that the curves remarkably 
show a nontrivial scenario described and reported by Ref. [3]. That is, in one direction 
in the sky a was smaller at the time of absorption, while in the opposite direction it 
was larger. More explicitly, we illustrate 3 pairs of Aa/cto(z)-predictions in table 1 
as examples. In the table, 0 = 0 (or n) means the quasar sight line is parallel (or 
anti-parallel) to direction of Q x , 'd\ D ^) = {17.4 h, —59°}. For each z, there are 
two values of (Act/ctojth with opposite signs, which just matches the expectations of 
observations. Such a theoretical picture is subtle. In addition, it were also reported as 
a dipole form in [3] 

Act - . 

-~ A obs cos 0, with A obs = (1.02 ± 0.21) x 10" 5 , (5.120) 

Ct 0 

where A 0 jj S means the observation value of amplitude A. Theoretically, figure [3 indi¬ 
cates that when z ~ 2 to 4, A t h ~ 1. x 10 -5 which coincides with A obs . 

4. We further plot the curves of Act/cto of (15.1141) with (15.1191) in more wide ^-region 
including radiation epoch of the Universe in Figured] (that epoch roughly corresponds 
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Figure 7: a-varying in the region of 0 < z < 3. 0 is angle between quasar sight line and axis Q 1 . 0=0 (or 
7 r) means the quasar sight line is parallel (or anti-parallel) to direction of <5 1 (^P^, = {17.4 h, —59°}. 

When z fixes, there are two values of Aa/ao(z) with opposite signs. Aa/ao(z) is given by Eg. (15.1141) 
with parameters Eg. (15.1191) . Three Keck’s data in 2004 reported and discussed by [9] 39] are shown for 
comparison. 


to z > 3 x 10 3 ). For 0 = 0, the radiation epoch limit of Aa/a is 4.7 x 10~ 5 , while for 
0 = 7r, this limit is about —5 x 10~ 6 . Therefore, we find that in that epoch the dipole 
form (15.1201) is no longer true. 


\a/a 0 



z(redshift) 


Figure 8: a-varying in the region of 0 < z < 4000. 


6 a -Varying in Whole Sky 


In the last section, the a-varying for 0 = {0, 7r} (or for case that both distant atom 
and quasar lie at Q 1 -axis) was studied and the model’s parameters R, M°, M 1 have been 
determined. Now we discuss general case of 0 G {0, 7r}, i.e., the case of Q 1 ^ 0, Q 2 ^ 0, 
and Q 3 ^ 0 (see Fig. [5]). The corresponding M-Beltrami metrics (I3.27P reads 



/ (q°—m°) 2 , i 
it 2 (<r( M )) 2 ^ TTW 

R2( ct (M) )2 

— (Q°—M°)Q 2 
R 2 ( ct ( m ) ) 2 

— (Q°—M°)Q 3 

\ R 2 ( CT ( M) ) 2 


R 2 ( ct (M ))2 

(Qi-M 1 ) 2 l 
R 2 { c t ( m )) 2 77 "! 

(Q 1 -M 1 )Q 2 
R 2 ( a ( M )) 2 
{Q 1 -M 1 )Q 3 
R 2 ( a ( M )) 2 


~(Q°-M°)Q 2 
R 2 ^ 1 )) 2 
(i Q 1 -M 1 )Q 2 
R 2 (^ M )) 2 

(Q 2 ) 2 _L_ 

i? 2( CT (M))2 a (M) 

Q 3 Q 2 

R 2 (a( M )) 2 


_(Q0_M°)Q 3 \ 

it 2 (cr("))2 

(Q 1 -M 1 )Q 3 
R 2 { c t ( m )) 2 

q 3 q 2 

R 2 ( rr ( M )) 2 
(Q 3 ) 2 l 

R 2 {c r ("))2 _ 57 mT / 

( 6 . 121 ) 
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where 


a (M) = x (3°(*) - My - (Q\z) - M 1 ) 2 - ( g 2 ( z )) 2 - ( g 3 (^)) 2 

= 1 - ~^[(Q°(z) - M 0 ) 2 - (Q(z) cos O - M 1 ) 2 - Q(z) 2 sin 2 Q\, (6.122) 

n. 

where Q 1 (z) = Q(z) cos© and ( Q 2 (z )) 2 + (Q 3 (z )) 2 = <5(z ) 2 — Q 1 (z) 2 = Q(z) 2 sin 2 @ have 
been used. From Eqs. fl6.12ip (16. 1221) . we have 

det£g°(g) = -(a (M) )- 5 . (6.123) 


From Eq. fl4.67p . we have 


Ai = 


1 + 


-^[(<3(z) cos© - M 1 ) 2 + Q(^) 2 sin 2 0]^ , A 2 = A 3 = -a (A/) . (6.124) 


We now focus on the derivations of a in this case. In section 5, we presented the proce¬ 
dure for calculating a step by step in detail based on b[!' 1 !\Q)\(q , z=qz=o)- Though the full 
b\„P ( Q) (16.1211) is more complex than B < ^J\Q)\^q2 = q3 =0 ) (15.781) . the calculations in section 5 
can be repeated smoothly. The resulting a z (Q)-expression (15.1051) keeps invariant except the 
{B$\Q), <x ( m) , Aj}| ( Q 2 = Q 3 =0) in the formula should be replaced by {B^\Q), a (A/) , A;}|(q^ 0 ) 
with (i = 1, 2, 3). Namely we have 


i 


a 2 (h2) = a 


V f t ( m ) 


+ 


(■ Q°-M °) 2 
R 2 c t( m ) 


■\J— det(H^ ) (Q))\/—A iA 2 A 3 ' 


with a = 


he 


(6.125) 


where a^ M \ det B^\Q), A* are given in Eqs. fl6.122p . (16. 123ft and (16. 124ft respectively. The 
corresponding ct-varying formula reads 


Aa 

«o 


_l_, (qQ-m°) 2 f 1 A 

R 2 C r( M ) \ i+vTl "!) 


iV 0 ^/-det(5g f) (Q))V-A 1 A 2 A 3 


- 1 , 


(6.126) 


where IVo is given in (15.110p . 

The a-varyings Aa/«o(^, ©) are shown in Figj9] by using (I6.126P in which the curves 
correspond to z from 0 to 4.5 and 0 = {0, 7 t/4 , 7 t/3 , 0.47T, 7 t/2 , 0.67 T 27t/3, 37t/4 , 7 t} 
respectively. We can see from the figure that: (i) When z were fixed, Aa/cco decreases along 
with 0 increases from 0 to 7r; (ii) In regions of {0 < 0 < 0.47 t} and {0.67T < 0 < 7r} ,a vary 
special spatially. That is, a could be smaller in one direction in the sky yet larger in the 
opposite direction at the time of absorption. This feature is consistent to the observations in 
Keck and VLT reported by [3] and ffl; («i) When 0 ~ {0.57T, 0.6(/r, 0 . 77 r}, the observation 
results of a-variations Aa/a 0 are nearly null. 

In order to show Aa/cco’s dipolar behavior more explicitly, we further plot figure of 
Ack/cko vs A r cos© in FigJTOl Three theoretical prediction curves corresponding to 0 = 
{0, 7t/4, 0.47t} and the experiment observation data reported by [3j are shown in the figure 
for comparison. It can been seen that the three curves are approximately close to each 
other in the region of r cos© = {—2.5GLyr, 2.5GLyr}, and their average gradient is about 
A ~ 1.0 x 10 _6 GLyr . This theoretical prediction value is consistent with the observation’s 
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Aa'/a ' 0 



z(redshift) 


Figure 9: a-varying Aa/ao(z,©) in the region of 0 < z < 4.5 and 0 < 0 < tt. The parameters 
{R, M°, M 1 } are shown in the figure. 



Figure 10: Curves of Aa/ao(z, 0) vs A rcosO. Two solidline curves and one dotted curve are shown. One 
of the solidline curves corresponds to O = 0 and the other is for O = 0.47T. The dotted line curve corresponds 
to O = 7t/4. The horizontal axis shows projection of atom’s “distance” r(z) = ct(z) onto Q 1 axis. t(z) and 
cosmology parameters (H 0: !!,«, Oa) were given in (12.161) and (12.171) respectively. The data reported by [3] 
are plotted with error bars. The parameters {R, M°, M 1 } are shown in the figure. 


(1.1 ± 0.25) x 10 6 GLyr reported in [3]. However, for absorbing systems with 0 ~ 0.47T and 
|rcos0| > 3.5GLyr, the curve with 0 = 0.47T in Fig JlTTI indicates that Aa/ao(0 0.4vr) ^ 
—Aa/ao(0 ~ (vr — 0.47t)). This means that the dipole term (i.e., the term cc cos0) in the 
expansion of Aa/ao is no longer dominating. For the absorbing systems with 0.47T < 0 < 
0.57T the situations is also similar. Observational experiments to check such predictions is 
called for. 


7 Summary and Discussions 

The spacetime variations of fine-structure constant a = e 2 /hc in cosmology is a new 
phenomenon beyond the standard model of physics. To reveal the meaning of such new 
physics is of utmost importance to a complete understanding of fundamental physics. The 
main conclusion of this paper is that the phenomenon of a-varying cosmologically with 
dipole mode dominating is due to the de Sitter (or anti de Sitter) spacetime symmetry in an 
extended special relativity called de Sitter invariant spacial relativity (dS-SR). Specifically, 
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the logic that leads to this conclusion are summarized as follows: 

1. The Keck+VLT data that imply varying a are results of measuring spectra of atoms 
(or ions) in distant absorption clouds. So it is legitimate to use the electron wave 
equation of atoms (typically, the Hydrogen atom) to discuss this issue. 

2. As usual, QM equations for spectra in atoms are defined in inertial coordinate systems 
to avoiding ambiguities caused by inertial forces. So, it is necessary to take the local 
inertial coordinate systems in FRW Universe for discussing both laboratory atoms and 
the distant. 

3. When Einstein’s cosmology constant A = 3/R 2 ^ 0, the metric in the local inertial 

coordinate systems in FRW Universe has to be Beltrami metric B^ u (x) or M-Beltrami 
metric but cannot be Minkowski metric r) lw . 

4. Since there exist both temporal and spatial variations for a in cosmology, B$?\x) is 
suitable. The de Sitter pseudo-radius parameter R and the Minkowski point parameters 
{M°, M 1 } are expected to be determined by fitting to the observations. 

5. As usual, dS-Dirac equation for hydrogen can always be reduced to spectrum equation 
of hydrogen. In this way [21J[22] both coefficient of Dirac-kinetic energy operator term 
— icl ■ V (i.e., h z (Q)c) and coefficient of Coulomb potential term — 1/r (i.e., e^(H)) are 
derived explicitly. Then a z (D) = e^(H)/ h z (Q)c and Aa/cto = (ay(H) — ao)/ao can be 
calculated. 

6. According to w, we focuss the best-fit direction about right ascension ~ 17.5 h, 
declination ~ —61 deg, and calculate Acc/cco in this region. Comparing the theoretical 
prediction with the observation results reported in [30], the model’s parameters are 
determined: R ~ 500GLyr, M° ~ —lOOGLyr and M 1 ~ —22GLyr. Surprisingly 
but not strangely, the amazing observational discover about Aa/a 0 -dipole in [30] is 
reproduced by this dS-SR theoretical model. This is a main result of this paper, which 
could be thought as a nontrivial evidence to support dS-SR. 

7. The a-varying in the whole sky have also been studied in this model with the same 
parameters. The results are generally in agreement with the estimations in [3j|4]. 

To get things more straight, let’s go back temporarily to the beginning of this paper again. 
If the theory of SR is exactly Einstein’s SR (E-SR) with metric r)^, the a will not vary over 
cosmic time, i.e., Acc/ct = 0. However, the observations among cosmological distances go 
directly against it, i.e., (A a/a) 0 b s ^ 0. Even though people would be open-minded to face 
this big challenge [2j, the simplest answer to the puzzle seems to be that E-SR might not 
be exact for cosmological spacetime scale, or E-SR needs to be extended. To the best of 
authors’ knowledge, the most natural consideration to extend E-SR in the framework of SR 
is the works due to Dirac(1935)-Inonii and Wigner(1953)-Gursey and T.D. Lee (1968)-Lu, 
Zuo and Guo (1974) [T4LIT8] . Those works produced the theory of de Sitter invariant special 
relativity (dS-SR). Using the dS-SR in this paper, we correctly work out the prediction of 
Acc/ccq which is consistent with the data reported by [210] . Hence the puzzle of a-varying 
over cosmological time could be considered solved. Our approach could be considered a very 
simple answer to the problem, if not the simplest answer. 
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The relativity principle problem in curved spacetime with constant curvature has been 
solved in dS-SR in na ng via introducing Beltrami metric (see [20] and the Appendix 
of [21]), and hence E-SR is the limit of dS-SR with \R\ —y oo. Since | it (-value determined 
in the present paper is fortunately cosmology huge (~ 500GLyr), we argue that such dS-SR 
would not contradict the experiments verified E-SR within the error bands. Furthermore, 
the many-multiplet (MM) method [5[[6] used in [|3[[4] is itself it -1 -free. That the it is so huge 
(or almost infinity) means that the estimates to A a/a in ^Jj4] are reliable approximately. 
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A Beltrami Metric and de Sitter Invariant Special Rel¬ 
ativity 

In this Appendix we briefly interpret the Beltrami metric and the de Sitter invariant 
Special Relativity (dS-SR). 

1. Beltrami metric: 

We derive the expression of Beltrami metric (II.9ft in the text. We consider a 4- 
dimensional pseudo-sphere (or hyperboloid) S A embedded in a 5-dimensional Minkowski 
spacetime with metric t)ab — diag(l , -1,-1,—1,-1): 

5 a : Vab£ A Z b = -R 2 , 

ds 2 = r] AB d^ A d^ B , (A. 127) 

where index A, B = {0,1, 2, 3, 5}, R 2 := 3A _1 and A is the cosmological constant. 5 a 
is also called de Sitter pseudo-spherical surface with radii R. Defining 

P 

:= -R-A-, with £ 5 ^ 0, and /i = {0,1, 2, 3}. (A.128) 

and treating P are Cartesian-type coordinates of a 4-dimensional spacetime with met¬ 
ric g IJiV (x) = B Ml/ (x), denoting this 4-dimensional spacetime as B A (call it Beltrami 
spacetime), we derive B^(x) by means of the geodesic projection of {5 a h- B a } (see 
Figure fill) . From the definition (1A.127I) . we have 

ds 2 = r]ABdi A di B \^A,B &Sh 

= ^d^dc - {d^f 

:= B^ u (x)dx^dx u . (A.129) 

Since £ A,B G 5a, and from (IA. 128j) and (IA. 127ft . it is easy to obtain: 

e = ^-e, ^ = hew +« 5 ) 2 = -^, 

R R a(x) 

£ £ 1 ri 'T»M/7'7* z/ 
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Figure 11: Sketch of the geodesic projection from de Sitter pseudo-spherical surface S\ to the 
Beltrami spacetime B\ via Ea. dA.128l) . 


where 


a(x) 


ri^x^x" 


Substituting them into Eq. flA.129ft . we have 


, 2 _ V^dx ll dx u (r] liu x^dx v ) 2 
a(x) R 2 a(x) 2 


B^(x)dx^dx u . 


Then, we obtain the Beltrami metric as follows 


Bfu/ (x) 


Vi*v Vnx x x r]^x p 

a(x) R 2 a(x) 2 


which is just Eq. (11.9ft in the text. 


(A.130) 


(A.131) 


2. Inertial reference coordinates and principle of relativity: 

The first Newtonian law is the foundation of the relativity. This law claims that the 
free particle moves with uniform velocity and along straight line. There exist systems 
of reference in which the first Newtonian motion law holds. Such reference systems 
are defined to be inertial. And the Newtonian motion law is always called the inertial 
moving law. If two reference systems move uniformly relative to each other, and if one 
of them is an inertial system, then clearly the other is also inertial. Experiment, e.g., 
the observations in the Galilco-boat which moves uniformly, shows that the so-called 
principle of relativity is valid. According to this principle all the law of nature are 
identical in all inertial systems of reference. 

Theorem 1\ The motion of particle with mass mo and described by the following 
Lagrangian 

L Newton = ^m 0 v 2 = -m 0 x 2 (A. 132) 

satisfy the first Newtonian motion law, or the motion is inertial. In (IA. 132ft . the 
Cartesian expression of the velocity is as follows 


v = x, and x = obi + x 2 j + x 3 k, 


(A.133) 






















where i-i = j- j = k k = 1 , and i j = i k = j k = 0. 

Proof: By means of the Euler-Lagrangian equation 

dL d dL dL d dL 

dx l dt dx l ’ <9x dt die 

(where d/dye = V := (d/dx 1 )i + (d/dx 2 )j + (d/dx 3 )k and etc) and L 
obtain 


(A.134) 


LNewton We 


x 1 = 0, x 1 = v l = constant , or x = v = constant. QED. (A.135) 


Theorem 2: The motion of particle in Minkowski spacetime described by 

y / rj tw dx> l dx h 


ds 

L E instein VTIqC 

dt 


-m 0 c- 


dt 


= — m 0 c 


(A.136) 


is inertial. 


The proof is the same as above, because both L Newton and L Einstein are coordinates x l - 
independent. Generally, any x-free and time t-free Lagrangian functions L(x) can 
always reach the result of (1A.135|) . However, when Lagrangian function is time- 
dependent that rule will become invalid. A useful example is as follows: 


L\(t, x, x) = —m 0 c : 


: V / 37A 1 


' 3/A(c 2 


x 2 x 2 + (x-k) 2 + c 2 (x-xt) s 


c 2 (3/A + x 2 — c 2 t 2 ) 2 


, (A.137) 


where a constant A ^ 0. The stick-to-itive readers can verify the following identity via 
straightforward calculations from (1A. 13711 : 


dL\ d dL\ 
<9x dt die 


x ■ 


d \ dL\ 


dye) die 


Noting that the Euler-Lagrange equation (IA. 134)1 reads 


dL\ d dL A d dL 


A 


d \ dL f 


x ■ 


dye dt die dt die 
and substituting (1A. 13811 to (1A. 139(1 . we have 

d \ dL A 


dye) die 


d \ dLf 


x • 


die J d± 


x 


T7T- =0. 


die J die 


Since 


d dL A 
die die 


= det 


d 2 L A 

dx l dxi 


7^0 


(A.138) 


(A.139) 


(A.140) 


(A.141) 


we have 

x = 0, x = v = constant , (A.142) 

which indicates that the particle motion described by Lagrangian function (IA. 13711 
is inertial, and the first Newton motion law holds. Thus, the corresponding inertial 
reference systems can be built. Noting 


hm L A L Einstein i 

A—>-0 


(A.143) 


it is essential and remarkable that a new kind of Special Relativity based on L A (I A. 1371) 
serving as an extension of the Einstein’s Special Relativity (E-SR) may exist. 
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3. de Sitter invariant Special Relativity (dS-SR): 

Following the Landau-Lifshitz formulation of Lagrangian |38] (see fl A. 13611 ). we examine 
the motion of free particle in the spacetime with Beltrami metric (1 A. 131 [) . From 
Eq. (12.18ft in text 


ds 

L dS — —m ° c yjy = 


(A. 144) 


we derive its expression in Cartesian coordinates. Setting up the time t = x°/c, B^ u (x) 
can be rewritten as follows 


ds 2 = B^ v {x)dx^dx 1 ' = good(ct) 2 + gij [( dx l + N l d{ct)){dx 2 + N 2 d(ct))] 

= c 2 (dt) 2 go 0 + g ij (^x l + N l ){^x J + N J ) , (A.145) 


where 


R 2 


fl'oo — 


a(x)(R 2 — c 2 t 2 ) ’ 

1 


— Vij l m 

9ij = 77 ^ + T 31~(„\2 VilVjmX x , 


N l = 


a(x) R 2 a(xy 
ctx 1 

R 2 - c 2 t 2 ‘ 


(A.146) 
(A.147) 
(A. 148) 


Substituting eqs. (lA.145l) - (lA.148l) into (1A.144I) . we obtain the Lagrangian for free par¬ 
ticle in B\. 

L d s = -m 0 c 2 \g 00 + + N^-xi + Ni). (A.149) 

V c c 


By using Cartesian notations (1A. 13311 and expressions of (1 A. 13011 (1A. 14611 (IA. 14711 
(1A. 14811 . the explicit expression of Lagrangian (1A. 149ft is: 


L d s = -m 0 c- 


R 4 


(R 2 + x 2 — c 2 t 2 )(R 2 — c 2 t 2 ) 

Z -/2 „ 2 + 2_2 


X 

~V + 


CCx" 


2t(x • x) 

R 2 + x 2 — c 2 t 2 c 2 ' (R 2 — c 2 t 2 ) 2 ' R 2 — c 2 t 2 


+ 


R 2 


(R 2 + x 2 — c 2 t 2 ) 2 \ c 


x • x ctx 2 
-+ 


R 2 - c 2 t 2 


1/2 


2 R 2 (c 2 — X 2 ) - X 2 X 2 + (x ■ ±) 2 + C 2 (x - xty 

= —rriQC R\ 


c 2 (R 2 + x 2 — c 2 t 2 y 


(A.150) 


where x 2 = (x ■ x). Noting R 2 = 3/A (see, e.g., Eq.(15) in Ref. [23]), and comparing 
L d s with L a (£,x, x) of (1A. 137ft , we hnd 


L d s = L\(t, x, x) = - m 0 c yj B^ v (x)x^x v , 


(A.151) 


which is the Lagrangian for free particle mechanics of dS-SR. Since (1 A. 1431) . when 
\R\ —* oo, the dS-SR goes back to E-SR. 
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4. de Sitter transformation to preserve Beltrami metric B^ u : 

In j2Qj (see Eqs. (41)-(43) in [20]), we have shown that under Lu-Zou-Guo (LZG) 
transformation (see also Eg. HB.1611) below) preserves Beltrami metric B^ u . When 
space rotations were neglected temporarily for simplify, the transformation both due 
to a Lorentz-like boost and a space-transition in the x 1 direction with parameters 
j3 = x 1 /c and a 1 respectively and due to a time transition with parameter a 0 can be 
explicitly written as follows: 


t —y t 

x 1 —y x 1 

x 2 —> x 2 
x 3 —>• x 3 



ct — fix 1 — a 0 + j3a 1 + 



x 1 — (3ct + /3a° — a 1 + 



a 0 —/3a 1 a 0 ct—a 1 x 1 — (a°) 2 +(a 1 ) 2 


R 2 



a 1 —/3a° a°ct—a 1 x 1 — (a°) 2 +(a 1 ) 2 



(A.152) 


where 7 = 1/ y/l — (3 2 . It is easy to check when R —y oo the above transformation goes 
back to Poincare transformation (or inhomogeneous Lorentz group ISO( 1, 3) transfor¬ 
mation) in E-SR. 


5. Conserved Noether charges of SO(4:, 1) of dS-SR: 

The external spacetime symmetry of dS-SR is SO( 4,1). According to Neother theorem, 
the corresponding 10-Noether charges are energy E, momentums p l , boost charges K l 
and angular-momentums L l . All have been derived in [20]. The results are as follows 


Noether charges for Lorentz boost : K l = moTc(x l — tx l ) 
Charges for space — transitions (momenta) : p l = m 0 Tx*, 

Charge for time — transition (energy) : E = m 0 c 2 T 

Charges for rotations in space (angularmomenta) : L l = e) k x j p k , 


where the Lorentz factor of dS-SR is: 


r 



1 


(x-x) 2 —X 2 X 2 


+ 


(x-xp 2 

R 2 


(A.153) 


(A.154) 


It can be checked that E = p l = K l = L l = 0 under the equation of motion x l = 0 (or 

x = 0) [2D]- 


B Modified Beltrami Metric and de Sitter Invariant 
Special Relativity 

We provide a brief introduction to Modified Beltrami metric (M-Beltrami metric) and 
the corresponding dS-SR. 

1. M-Beltrami metric: Eqs. (11 . 10p and (11.lip are the definition of M-Beltrami metric 
B$?\x). Being different from B liu (x), the coordinate components of Minkowski point 
for b\^!\x) is instead of the origin of spacetime system xC Introducing notation 

y^ = x^- AP (B.155) 
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then 


B\!Z\x) = B„„M- 


(B.156) 


The Landau-Lifshitz action is 


S = —me J \JB l ' p A J\x)dxf 1 dx l/ = —me J \JB pu (y)dy^dy u , (B.157) 

where dM p = 0 were used duo to constancy of MR The Lagrangian Lm-cls reads 

Lm-<is = —mc\J B^\x)x^x u . (B.158) 


Then, from 5S = 0 and Eq. (lA.142|) . we have 

<% 0 ) 2 ’ 

where y = x — M, y° = ct y = ct — M° (see Eq. (|B.155l) ). Eq. (IB.159j) becomes 

d 2 x 


(B.159) 


dt 2 


= x = 0, 


(B.160) 


which means that the free particle moves with uniform velocity and along straight line 
in the dS-SR based M-Beltrami metric. Consequently, the first Newtonian law holds 
for Lm-<is Eq. (IB. 1581) and inertial coordinate frames are well defined. 


2. Spacetime symmetries of M-Beltrami metric and the motion integrals. In [20] (see Eqs. 
(41)-(43) in [20]), we have shown that under Lu-Zou-Guo (LZG) transformation 


r 


LZG ) ~fi 

= ±cr(a) 1//2 cr(a,:r) 

D t 

= L p + R^ 2 p up a p a x 

L : 

= (L p )eSO( 1,3), 

o{y) 

= 1 - , 

a(a,y) 

= 1 - -^rj pu a p y u , 


(B.161) 


the Beltrami metric transformation reads: 


B,Jy) 


LZG 


> B , 


dy x dy f 


fJ,V \ 


dyv dy 

(IB.1621) leads to the invariance of action of (jB. 157j) : 

s -^^4 s — s. 


B\ P {y) = B : 


LLV v 


(B.162) 


(B.163) 


The corresponding Noether chargers or conserved motion integrals are as follows: 


Noether charges for Lorentz boost : K l = mTc(y l — t y ^~) 
Charges for space — transitions (momenta) : p l = mT^-, 

Charge for time — transition (energy) : E = mc 2 Y 

Charges for rotations in space (angularmomenta) : L l = d- k y 3 p k , 


(B.164) 


32 

























where the Lorentz factor of dS-SR is: 


1 


r = 


1 C 2 ( dty ) + C 2 i ? 2 


. d y y 

dty I 


dy V 

dty J 


(B.165) 


+ J_ ( v - t 

+ R 2 [ y l V dty 


Using Eq. flB.155p . we have the expressions in x frame: 


K 1 = mTc[x { — M i — (t — M°/c)x% (B.166) 

p l = mTx 1 , (B.167) 

E = me 2 r, (B.168) 

U = e) k (x-My P \ (B.169) 


and 


r 


1 



[(x—M)-x ] 2 —(x—M) 2 x 2 

^PR 2 


+ 


[x— M— -k.it— M°/c)] 2 

W- 


(B.170) 


It is straightforward to check that E = p 1 = K l = L l = 0 under the equation of motion 
x l = 0 (or x = 0) and M 11 = const.. 
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